Gravitational Waves in the Spectral Action of Noncommutative Geometry 
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The spectral triple approach to noncommutative geometry allows one to develop the entire stan- 
dard model (and supersymmetric extensions) of particle physics from a purely geometry stand point 
and thus treats both gravity and particle physics on the same footing. The bosonic sector of the 
theory contains a modification to Einstein-Hilbert gravity, involving a nonconformal coupling of 
curvature to the Higgs field and conformal Weyl term (in addition to a nondynamical topological 
term). In this paper we derive the weak field limit of this gravitational theory and show that the 
production and dynamics of gravitational waves are significantly altered. In particular, we show 
that the graviton contains a massive mode that alters the energy lost to gravitational radiation, in 
systems with evolving quadrupole moment. We explicitly calculate the general solution and apply 
it to systems with periodically varying quadrupole moments, focusing in particular on the the well 
know energy loss formula for circular binaries. 
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I. INTRODUCTION 

Noncommutative Geometry (NCG) is a gravitational 
theory which, even in its simplest form, can explain the 
Standard Model of particle physics, and account for all 
current experimental data, in a rather simple and cer- 
tainly elegant way. The simple — in the sense that it 
generalizes the continuum Riemaniann manifold by con- 
sidering its product by a discrete two points space — 
NCG proposal should be certainly replaced by a less triv- 
ial noncommutative space as one reaches Planckian en- 
ergy scales. Nevertheless, this is the proposal we have at 
hand, and given its success in accounting for the Stan- 
dard Model of particle physics, it offers a valid theoretical 
framework to address early universe issues. Alternatively, 
one can use experimental high energy physics data and 
astrophysical observations/measurements in order to test 
this NCG proposal and constrain its parameters. This is 
the approach used in this study. 

One should indeed view this NCG proposal as an ef- 
fective theory, which can however offer a valuable infor- 
mation about any NCG approach. In what follows we 
derive explicitly the weak field limit of this gravitational 
theory and then show that the production and dynam- 
ics of gravitational waves are both considerably modified 
from those obtained within the familiar General Relativ- 
ity approach. 

More precisely, in Section |lT] we give a short introduc- 
tion to the noncommutative geometry spectral action, 
the framework within which we will then focus. In Sec- 
tion |lll] we first state in detail the conventions and sig- 
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nature we use and we analyze the issue of gauge con- 
ditions. We then analyze linear perturbations around a 
Minkowski background metric and we solve the noncom- 
mutative geometry gravitational wave equation in terms 
of the retarded Green's function. We find that gravita- 
tional waves are only sourced from systems with a non- 
trivial quadrupole moment, as within General Relativity, 
while the NCG theory contains massive as well as mass- 
less gravitons. In Section ITVl we concentrate first on some 
simple and then on some physical examples. Using the 
requirement that the mass of the gravitons must be posi- 
tive (and real) , we can fix the sign of the couplings in the 
NCG spectral action approach. We then calculate the 
energy loss for a circular binary system and compare it 
to the results obtained from standard General Relativity. 
We conclude that the amplitude of modifications within 
NCG is small, nevertheless the NCG approach leads to 
some distinctive features which we analyze. We round 
up with our conclusions in Section Fvl 



II. NONCOMMUTATIVE GEOMETRY 
SPECTRAL ACTION 

In the NonCommutative Geometry [l|, [2| approach, the 
Standard Model (SM) of electroweak and strong interac- 
tions is considered as a phenomenological model, which 
dictates the geometry of space-time, so that the asso- 
ciated Maxwell-Dirac action functional produces the SM 
with all known experimental results. The outcome of this 
approach is a geometric space defined by the product, 
jM X J^, of a continuum compact Riemaniann manifold, 
AA, and a tiny discrete finite noncommutative space, J-, 
composed of only two points. Such an almost commuta- 
tive space is the simplest extension of the more familiar 
commutative space upon which General Relativity is for- 
mulated. Certainly one should not expect the validity 



of this simplistic approach to hold at the Planck scale, 
which is the scale at which all notion of classical geometry 
loses its meaning. 

The metric dimension of the product geometry Ai x 
J^ is 4, the same as the ordinary space-time manifold. 
Thus, the metric dimension of the noncommutative space 
T is zero, while for noncommutative spaces one must 
distinguish between the metric dimension and the KO- 
dimension. The internal space J^ has XO-dimension 6 
to allow fermions to be simultaneously Weyl and chiral, 
whilst it is discrete to avoid the infinite tower of massive 
particles that are produced in string theory. 

The noncommutative nature of T is given by the real 
spectral triple {A, H, D) that generalizes Riemannian ge- 
ometry to the noncommutative setting; A is an involution 
of operators on the finite-dimensional Hilbert space % of 
Euclidean fermions, and D is a, self-adjoint unbounded 
operator in %. The choice of Hilbert space has no im- 
portance, since all separable infinite-dimensional Hilbert 
spaces are isomorphic. The algebra A, related to the 
gauge group of local gauge transformations, is the alge- 
bra of coordinates. A space is described by the algebra of 
coordinates, which in the context of NCG is represented 
as operators on a Hilbert space. Since real coordinates 
are represented by self-adjoint operators, all information 
about a space within NCG is encoded in the algebra of 
coordinates A. By assuming that the algebra constructed 
in A^ X J^ is symplectic-unitary, A must be of the form 



A^Ma{M)®Mk{C) 



(1) 



k = 2a, H is the algebra of quaternions. The choice 
fc = 4 is the first value that produces the correct number 
of fermions in each generatioii, i.e., fc^ = 16 fermions in 
each of the three generations |3J . 

The operator D corresponds to the inverse of the 
Euclidean propagator of fermions, and is given by the 
Yukawa coupling matrix which encodes the masses of the 
elementary fermions and the Kobayashi-Maskawa mixing 
parameters. The commutator [D,a\, with a € A, plays 
the role of the differential quotient da/ds, with ds the 
unit of length. The familiar geodesic formula 



d{x, y) = inf / ds 



(2) 



where the infimum is taken over all possible paths con- 
necting X to 2/, which is used to determine the distance 
d{x, y) between two points x and y within Riemannian 
geometry, is replaced by 

d{x,y) = sup{|/(a;) - f{y)\ : / e AI|[A/]|| < 1} , (3) 

where D is the inverse of the line element ds, within the 
noncommutative spectral geometry. 

The fermions of the SM provide the Hilbert space Ti, 
of a spectral triple for the algebra A, while the bosons of 
the SM, including the Higgs boson, are obtained through 
inner fiuctuations of the Dirac operator of the product 



M. X F geometry. Hence, the Higgs boson, which gen- 
erates the masses of elementary particles through spon- 
taneous symmetry breaking, becomes just a gauge field 
corresponding to a finite difference. Note that the cor- 
responding mass scale specifies the inverse size of the 
discrete geometry J- . 

Applying the spectral action principle, according to 
which the action functional on spectral triples depends 
only on the spectrum of the line element, i.e., the inverse 
of the Dirac operator, to the inner fiuctuations of the 
product geometry Al x J", one recovers the SM coupled 
to gravity in the Euclidean form. Thus, the NCG spec- 
tral action approach — limited to the classical level even 
though it can a priori be quantized — offers an elegant 
geometric interpretation of the SM , the most successful 
phenomenological model of particle physics. 

To be more precise, the SM Lagrangian — including 
mixing and Majorana mass terms for neutrinos, mini- 
mally coupled to gravity — can be successfully recov- 
ered from the asymptotic expansion of the spectral action 
functional 



H'ij)) ■ 



(4) 



where / is a positive even function of the real variable 
and A fixes the energy scale. Note that D/A is dimen- 
sionless since the Dirac operator, being a differential op- 
erator, has dimensions of mass. The physical Lagrangian 
is thus obtained from the asymptotic expansion in the 
energy scale A of the spectral action functional, Eq. (j4]). 
More precisely, using heat kernel methods one can write 
the square of the Dirac operator in terms of the inverse 
metric, the unit matrix and two matrix functions com- 
puted from D and show that the trace, Eq. (j4|) above, 
can be expanded in a power series as a function of the 
inverse scale A and it can thus be written in terms of the 
geometrical Seeley-deWitt coefficients a„, as [ll[ 



E ^4^nA4 



(5) 



n=0 



where the function F is defined such that F{D'^) — f{D). 
Defining the moments 



fk 



f{u)u''-'du , for /fc > , 



(6) 



and /o = /(O), one finds 

Fi = 2/4 
F2 = 2/2 
Fq = /o 



F„ 
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for 71 > 1 



(7) 



while the coefficients a„ are known for any second order 
elliptic differential operator. 



The coupling with fermions can be obtained by includ- 
ing an additional fermionic term 



1 



< J^p,Dip > 



(8) 



in Eq. ^, where J is the real structure on the spectral 
triple and ip is a spinor in the Hilbert space H of the 
quarks and leptons. 

The spectral action approach leads naturally to the 
merging of the three coupling constants at the unification 
scale, 52 = ff3 = y/5/3gi , it provides neutrino masses and 
mixing as well as the see-saw mechanism, and it predicts 
a heavy Higgs mass. 

The spectral action, Eq. ([4]), can be expanded in pow- 
ers of the scale A in the form 



J2 fkA''-f\D\-''+fmD{0)+O{l) , 



(9) 
where fk are the momenta of the function / given in 
Eq. ([6|), the noncommutative integration is defined in 
terms of residues of zeta functions, and the sum is over 
points in the dimension spectrum of the spectral triple. 

The physical Lagrangian that one obtains in this ap- 
proach, contains, in addition to the full SM Lagrangian, 
the Einstein-Hilbert action with a cosmological term, a 
topological term related to the Euler characteristic of 
the space-time manifold, a conformal Weyl term and a 
conformal coupling of the Higgs field to gravity. Note 
that the coefficients of the gravitational terms depend on 
the Yukawa parameters of the particle physics content. 
Within the NCG spectral action, one works in Euclidean 
rather than Lorentzian signature, assuming that one can 
get back to the Minkowski signature through Wick rota- 
tion. 

One then sets the parameters of the NCG spec- 
tral action at the (unification) scale A; predictions at 
lower energies are recovered by running the parameters 
down through Renormalization Group Equations (RGE) . 
Hence, the spectral action at the unification scale A of- 
fers a framework to investigate early universe cosmologi- 
cal models [6l-ll(j|. while extrapolations to lower energies 
can be obtained via, firstly, RGE and secondly, inclusion 
of nonperturbative effects in the spectral action. 

Adopting Euclidean signature, the gravitational part 
of the asymptotic formula for the bosonic sector of the 
NCG spectral action, including the coupling between the 
Higgs field (j) and the Ricci curvature scalar i?, is [J| 
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Note that H is a rescaling H = (-\/a/o/7r)0 of the Higgs 
field (j) to normalize the kinetic energy; the momentum 
/o is physically related to the coupling constants at unifi- 
cation and the coefficient a is related to the fermion and 



lepton masses and lepton mixing. 

In the above action, Eq. (|10p . the first two terms only 
depend upon the Riemann curvature tensor; the first is 
the Einstein-Hilbert term with the second one being the 
Weyl curvature term. The third term 



R*R* = \e 



^a/3^5 



jyafi rj^S 
^^^u -'^pa- 



is the topological term that integrates to the Euler char- 
acteristic and hence is nondynamical. The spectral action 
contains one more term that couples gravity with the SM, 
namely the last term in Eq. ([TU|) . which should always 
be present when one considers gravity coupled to scalar 
fields. This coupling can have significant consequences 
at high energies, such as in the early universe |6l4lC|. 
however in this paper we will be concerned with the low 
energy weak curvature regime where this term is small. 

Neglecting the nonminimal coupling between the Higgs 
field and the Ricci curvature, the equations of motion 
derived from the Lorentzian version of spectral action 
above read fal 



R^"" - \9^''R - 327rGao 



2^.1^;" + C^^'"'i?A« 



= STrGT^attcr : (11) 



implying that the NCG corrections vanish [0] for 
Friedmann-Lemaitre-Robertson- Walker (FLRW) cos- 
mologies. [The reader is directed to subsection IHII A for 
a definition and discussion of the Lorentzian conventions 
used.] 

We will be concerned with linear perturbations around 
a Minkowski background metric in the synchronous 
gauge, so that the perturbed metric reads 

.9^. - diag ({a(i)}2 [-1, (5,.^ + h,^ (x))]) , (12) 

where a{t) is the cosmological scale factor. Throughout 
this paper we work in a flat background and hence a(t) — 
1 and d = dt/dt = 0. The remaining gauge freedom can 
be completely fixed by setting Vi/i*-' = 0. [A detailed 
discussion of this gauge fixing is given in subsection IHIC 
below.] 

In Section [HT] below we show that, the linearized equa- 
tions of motion derived from the NCG spectral action, 
for such perturbations, read 



(D - /3^) Dh'"' 



o2 ^^T^G 

P „/l mat- 



(13) 



where T^attor i^ taken to lowest order in h^"^ . This implies 
that it is independent of h'^'^ and satisfies the conserva- 
tion equations 



d 



T^^=0 



(14) 



It is important to note that /?, defined as 

/32 = -l/(32^Gao) , (15) 



in Eq. ()13p . turns out to play the role of a mass and hence 
has to be real and positive, implying that ao < 0. In the 
following we will see that, for ao > 0, the gravitational 
waves evolve according to a Klein-Gordon like equation 
with a tachyonic mass, and hence the background, which 
in our case is Minkowski space, is unstable. We can thus 
restrict to ao < for Minkowski space to be a (stable) 
vacuum of the theory. 



III. PERTURBATION EQUATIONS 

To write down the linearized equations of motion, we 
will first discuss our conventions for the metric signature 
and the Ricci tensor. 



A. Conventions for 7?^,^ and Signature 

In this paper, we are using conventions in which the 
signature is (—,+,-(-, +) and the Ricci tensor is defined as 
Rp.i' = R''p.^p, with Rf^^p'^ujcr = [Vp, Vi'l^p- In General 
Relativity such choices are merely conventions, which are 
relatively unimportant (provided of course that one is 
consistent), here however the situation is very different. 
The Lorentzian version of the NCG action, Eq. (fTO|) . that 
we use reads 



1 



i? + «oC.,„,C^''''"+Toi?*i?* 
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- ^oi?|Hp) V^d^a; . (16) 

It is thus clear that the conventions used to define, for 
example, Rfj,i, will radically alter the theory, unless one 
also alters the (signs) of the couplings. Specifically, con- 
sider using the opposite convention for the Ricci tensor, 
which introduces a negative sign on all terms depending 
on Rfi^, but not on terms depending on R^^pa- Since 
our action now contains terms of both kinds (i.e., R and 
Cfj.i^paC^"'"^), this change of convention introduces a rel- 
ative sign change. This can simply be compensated for 
by changing ao — >■ —ao- However, without this change 
the action is very different. An exactly similar change 
happens if we considered a different choice of conven- 
tion for the signature or the sign of R^^pa, as these 
both introduce a sign change in the R term, but not 
the C^^paCf'^P" term. 

As shown in Refs. |a-l^, the presence of the non- 
minimal coupling of curvature to the Higgs field can have 
significant effects of the cosmological dynamics and one 
may wonder whether the sign ambiguities discussed here 
may affect these results. Fortunately, from Eq. (IT51) it is 
clear that the relative sign between the Einstein-Hilbert 
term (i?) and the nonminimal coupling (i?|Hp) is inde- 



pendent of any convention (since they both contain R). 
In a cosmological setting (i.e., for FLRW geometries) the 
Weyl term vanishes and hence, in homogeneous cosmolo- 
gies, the only NCG affects come from the nonminimal 
coupling (^], allowing such issues to be avoided. 

Since the underlying NCG theory is only developed for 
the Euclidean signature, it does not provide a guide for 
the Wick rotation to the Lorentzian space. Hence the 
choice of the sign of the couplings, appropriate for a par- 
ticular choice of convention, can only be made by testing 
the physical consequences of the theory. In the follow- 
ing section we will show that gravitational waves offer an 
excellent probe of the couplings in this theory, but even 
without an in-depth analysis, from Eq. (J13p one can im- 
mediately see that the coupling ao must be negative (in 
the conventions used here). If it were not, the (3 param- 
eter would be complex and this would correspond to a 
tachyonic mode of the graviton (we refer the reader to a 
discussion below). This would indicate that Minkowski 
space-time is unstable to small perturbations. If we rule 
out such a possibility on physical grounds (or require 
that the Lorentzian version of the NCG action should 
admit Minkowski as a stable vacuum) then we can re- 
strict ao < with the conventions used here. Consider 
for example the consequence of changing the signature, 
so that D — >■ — n. Such a change of convention essen- 
tially changes the sign of 0^ and the conclusions would 
be reversed^. 



B. Linearized Equations of Motion 

Variation of the gravitational part action Eq. ([2 
w.i.t. the metric gp^y leads to the following addition to 
the Einstein tensor of General Relativity (GR), 



^NCG 



n Einstein 

+2ao (2 Va VkC^""^ + C^^^^R^, 



(17) 



where as usual 

Given the convention used here to define the Ricci tensor 
the Weyl Tensor is explicitly given as, 

CpXuK = Rtj.XvK + {9ii[vRk\\ — 9\[uRk.\ij) (18) 
-l;9^l.[ugK\\R■ 



^ The concerned reader should note that exactly the same situation 
arises in standard Klein-Gordon equation, where the sign of the 
mass term is changed under a change of signature. 



Using the contracted Bianchi identity 

V^^mAi^k = -(VAi?pi.- Vm-Rai.), (19) 

and its remaining trace 



V'' ^ak — ^ Va R, 
we can arrive at the following expression 



(20) 



2 V V'^C'a'Ai/k — —Cxf^K^R '^ — V Va {R^h' — t9p.i^R) 

1 2 

+ o Vp Vi>^ ^ "^RiipR^v + ■::RRvti, 

+\9AR-^R^'' -\r^) ■ (21) 

Notice that the expression above for 2 Va V^C'^''""^ shows 
that the C^'^^'^Rf^x term in Eq. ^} exactly cancels in 
favor of terms which are of second order of solely the 
Ricci tensor and/or Ricci scalar. 



We now follow the standard procedure of perturbing 
about a flat metric, where 

9,.u = ^l^..+l^.. , 5^"^ = t?^"^ - 7^^ (22) 



and 



7 = 7% = ^^^^7^. 



(23) 



all tensor indices are raised and lowered using the back- 
ground metric 77^^ (except for the indices of 5^^ and g^'^). 
To first order in metric perturbations we then have 

2{vaV«C^^"'^^} = 

-i (ry^'^.a" - d''d'')dxd^^ + 0(7^) . (24) 

where 0(7^) denotes any second order combinations of 
7^1/ and we have defined 



Itiu = Ip.^ - 2^M''7 , 



(25) 



i.e., the trace reverse of 7^^. 



Similarly, to linear order in metric perturbations the 
Einstein tensor is simply 



G^rnstein = +^aAaV''+0(7') 



(26) 



C. Gauge Conditions 

In calculating [6| the linearized equations of motion, 
the traceless transverse gauge was imposed on the metric 
perturbations hfi^] here we explicitly show that this is 
indeed a valid choice. As before, we denote metric per- 
turbations that have not been gauge fixed by 7^^ and 
reserve h^j^i, for the final, gauge fixed perturbations that 
correspond to the physical gravitational waves. 

As always we have a freedom due to diffeomorphism in- 
variance of the action to restrict the gauge of the metric 
perturbations. Explicitly, under a diffeomorphism gener- 
ated by ^^ the metric perturbations 7^1/ transform as 



old ^ 

! pU 



7;r = 7;l'+9^e.+9.e^. (27) 



Without loss of generality, one can impose the Lorentz 
gauge conditions 



d^i, 



flu 



0, 



(28) 



restricting the perturbations to be transverse, where we 
introduced the "trace reverse" of 7^1/, as in Eq. ([25]) . 

Choosing this gauge (and dropping the label new on 
7^1,), Eq. ([24]) simplifies to 

2{^xV.C^^'"'}^ 

-^d^d^^i^Oxd^r" + liv^^^dxd^ - d''d'')j) . (29) 

Combining the above result with the Einstein contribu- 
tion, Gg^jjgjj,jjj, to the equation of motion, the left-hand 
side of Eq. (ITTI) is, to first order in 7^^, given by 

^NCG = ^'^Einstcin ^ 1^ \^ V k VaC''" '"' j 

+ ^d^d''(^dxd^r'' + liv^^^dxd^ ~ d'^d'^)^) . (30) 

However, the Lorentz gauge does not uniquely fix all the 
gauge freedom. More precisely, we still are free to per- 
form gauge transformations, generated by any ^^, that 
satisfy 



dad''^, = 



(31) 



since this still preserves the gauge condition, Eq. ([28]) . as 
can be checked directly. 

We can use this transformation, Eq. (|3T|) . to set (in 
the new frame) 7 = and 7oi = (i = 1, 2, 3) by solving 
the corresponding equations for ^i and their time deriva- 
tives on some initial surface t = to where no sources are 
present and further extending into a source free region 
(i.e., Tf^'^ = 0) [i2]- After performing these gauge trans- 
formations the source free equations of motion (i.e., the 



left-hand side of Eq. pTj) ) read 

- \^.^'^i^.. + ^d^d-dxd^-i^, = , (32) 
where we have made repeated use of the fact that 

Up to this point, the gauge restrictions 7 = and 
7oi = {i — 1, 2, 3) are the same as those typically used 
in General Relativity, which is to be expected since all we 
have used is the diffeomorphism invariance of the action. 
However, in determining whether we can set 700 = 0, 
the equations of motion are used and hence one might 
expect that this gauge condition will be different than 
that of General Relativity. To confirm its validity note 
that Eq. (UHl) implies 



57 



00 



dt 



0. 



(33) 



Now, using the equation of motion in the presence of 
matter source we arrive at 



V^ 700 - -gj VM V^ 700) 



1-^v^ 



V^7oo 



2k. 
2k 



(34) 



Recall that General Relativity is recovered in this set- 
ting by taking /3 — >■ 00. Thus, one can see that in this 
limit the equation simplifies to v 7oo = — (167rG')/c^roo, 
which fixes 700 to be a constant (assuming the space- 
time is asymptotically flat) away from the source. Fi- 
nally, a redefinition (gauge transformation) allows us to 
set 700 = 0. 



From Eq. (p4)) we see that, away from the source, 
V 7oo = is still a solution and hence we can fix 700 = 0, 
however this is no longer the only solution to Eq. (p4| . 
In particular, away from the source one could fix 700 via, 



1 2\ 

1-^ V J7oo = 0, 



(35) 



which clearly solves Eq. (|34)) . This would result in a 
modification of what is often referred to as the radiation 
gauge. In the following, we choose 700 = so as to be 
able to directly compare our results to the standard ones 
obtained within General Relativity. 

In cases where sources are present, the NCG equation 
of motion with gravity and normal matter is 






(36) 



condition of the radiation gauge ^ . But the explicit pres- 
ence of the trace 7 in Eq. p6p above can be eliminated 
by formally defining the tensor hfj,„ as, 

hf,^='ff,r,- -—^0^^{r]f,r.D~df,d^)j, (37) 

where the operator O is given by 



0= 1 - 






(38) 



This is a modification of the trace reverse of the metric 
perturbations that is usually used, however it performs 
the same task, namely removing the trace from the equa- 
tions of motion, Eq. (|36|| . As it can be easily checked, 
the Lorentz gauge condition, Eq. (pS)) . is satisfied by hfii, 
as long as it is also satisfied by 7^,^. Note that the trace 



of h^„ is 



^ = -'1 + ^^)7 



(39) 



so we see that indeed, this reproduces the trace reverse 
of "f^i/ in the limit /3 — >■ 00. Clearly then, when we are 
away from a source, we can impose that 7 = and this 
implies that h = 0. 

In terms of h/^^i, the equation of motion, Eq. p6p . is 



('-» 



Dh'"' 



167rG 



rpfiiy 



(40) 



/32 / c- 

Dropping the over-bars, this is exactly Eq. p^ . 

D. Green's Function 

The general physical solution to Eq. PO]) is given by 
h^"" = 2/3^^ / dSix')Gn{x, x')T^"'{x') , (41) 



where the Green's function Gr(x, x') satisfies the fourth- 
order partial differential equation: 

(D-l3^)DGn{x,x') = AirS^'^^ix - x') , (42) 

where the operators D above are acting on x. In order to 
find a solution Gb.{x,x') to Eq. (|^^ consider two distri- 
butions gi and 52 which satisfy the following second-order 
equations: 



(n-/32)gi = 471(5(4) (x - x') 
ng2 = AttS^'^^x - x') 



(43) 
(44) 



Since T^^, 7^ we are not free to impose the traceless 



^ Nor would this be possible if one had chosen the modified radi- 
ation gauge implied by Eq. I I35I I. 



Then one can easily verify that the combination 



1 



Gii{x,x') = -^{91-92) , 



(45) 



wih be a solution to Eq. (|42p . Physically we are interested 
in the retarded Green's function solution to Eq. (H^ . that 
is of the form 



Guix, x) = Q{t — t')g{x — x) , 



(46) 



where Q{z) is the Heavyside step function. So 91 and 
92 in Eq. (J45|) above must each be retarded solutions 
to Eqs. (j43|) and (J44|l . respectively. The other three 
combinations of retarded and advanced solutions violate 
causality and will not be considered further. So we have 
reduced the problem of finding the retarded Green's func- 
tion solution of the fourth-order, Eq. (|42l) . to finding the 
retarded Green's function solutions of the two second- 
order differential equations, Eqs. (|43)) and p4)) . The ex- 
plicit calculations are given in the Appendix and result 



ffiR = ^^^^lMe(cr-|R|)-2<5(r2)j>e(T),(47) 



92K - -25{t')Q{T), 

where we have defined 

T = t~t' 
R = r r' 



(48) 



r ^ ^{cTf - |R|2 , 

and Ji{x) is the first order Bessel function of the first 
kind. Thus, Eq. (|45]) implies that 



G^{x,x') = ^ZllMe(cT-|R|)e(T). (49) 

Note the absence of the delta function singularities on the 
light cone, consistent with the general analysis detailed 
in Ref. [20|. 

Using Eq. (|41|) one finds that the field is given by 



{cTf - |RP 



x^i(^/3V(cT)^-|RPJT^'^(r',t') 

xe(cr- |R|) . (50) 

One thus finds that the field is sourced only by regions 
within our past light code (i.e., cT > |R| ), which is 
expected for the propagation of a (positive-real) massive 
field. Notice that if /3^ < 0, corresponding to ao > 0, we 
find that an observed field is sourced from regions with 
space-like separation. This is due to the fact that /3^ < 
corresponds to a tachyon (complex mass) mode of the 



gravitational wave. 

If we consider the far-field limit, i.e. |r| w |r — r'|, we 
can write the spatial components of this field as 



h"' (r, t) 



2Gf3 



t--M 



dt' 



c^ {t - t'f - \v\^ 



X Ji (B^cm~t'f-\vA &'' [t') , (51) 

where we have, as usual, introduced the quadrupole mo- 
ment, 

D''=(i) = 4 /"dra;'x*'T"0(r,i) , 

the second time derivative of which is given by 



(52) 



drr'=(r,i) 



1 d^ 
6dt2 



D'^'it) 



(53) 



In conclusion, just as in the case of General Relativ- 
ity, we find that the gravitational waves are only sourced 
from systems with a nontrivial quadrupole moment. This 
is essentially due to the conservation of the energy mo- 
mentum tensor, which is unaltered in this theory. How- 
ever, the propagation of gravitational waves is signif- 
icantly altered, by the presence of additional massive 
modes. 



IV. EXAMPLES 

As a simple pedagogical example let us consider the 
case of I?*'^ « constant, for which one can explicitly per- 
form the integral in Eq. (j5ip to find 



W-(v,t) 



D 



■ik 



2G 

3c2|r| fixed 

X [sinh(^|r|) - cosh(^|r|) + 1] . (54) 



Thus, one recovers the standard result in the /3 — >■ 00 
limit, namely 



lim h^^{v,t) 

/3— ^00 



lim 



2G 



D 



ik 



(1 

fixed V 



-/3|r| 



3c4|r, 
(G^)/i"=(r,t) lim fl-e-'^l-^l), (55) 



where (GiR)/j»fe denotes the field in the General Relativis- 
tic case. This is expected, since the /3 — > 00 limit corre- 
sponds to taking ag — ^ and, as it can be seen from ei- 
ther Eq. (flO)) or Eq. (TTSt , one then recovers the Einstein- 
Hilbert result. 

The next simplest example is a system with a period- 
ically varying quadrupole moment, i.e. 



D'^ (t) ^ A'^ cos {uj'H + (1)'^) 



(56) 



where A^^ is a constant in time, w*-' is the frequency of 
the oscillations of the ij component and </>'■' is the phase, 
both of which we consider to be time independent. No 
summation over i, j is implied. Using this in Eq. (151 
one finds 



h'^ = 



AG(3A'3uj^- 
3^4 



sm{u:'H + f^)f4p\vl — 



+ cos {uj''H + (j>'^)f J /3\r\, 



/3c 



(57) 



where again no summation is implied and we have defined 
the functions. 
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III 



/s (2:, Z) 

/c (a;, z) 



y/s^ + x 



J\ (s)sin izys^ + x'^\ 

, „ „ Ji (s) cos [z\J s^ +; 
V s + x'^ ^ 



(58) 



These functions are highly oscillatory, with somewhat dif- 
ferent behavior for z > \ and z < 1. Because they have 
a typical frequency of the order of z, which is (in gen- 
eral) different to uJ^^ , the wave- form of the gravitational 
radiation, Eq. (|51l) and its time derivative, Eq. ([57)1 can 
experience beat phenomena. In particular, interference 
between the various functions can result in a significant 
enhancement of the amplitude. 

As a specific example, of significant physical interest, 
consider a pair of masses m\ and rui^ in a circular binary 
system, under the assumption that the internal structure 
of the bodies can be neglected. For such a system, or- 
biting in the xy-plane, one finds that the only nonzero 
components of the quadrupole are [l3| , 

b^^{t) = 12M|p|2sin(2V'(t))w3 



= -byy{i) 



D'^y (t) 

JJZZ 



-12^|ppcos(2V'(t))w^ 



(59) 



where /i = mim2/{mi + 7712) is the reduced mass of the 
system, \p\ is the magnitude of the separation vector be- 
tween the bodies, which is constant for circular orbits, ip 
is the angle of the bodies relative to the x-axis and u! = 4' 
is the orbital frequency, which for this simple system is a 
constant given by 



c^ = 1^ == IpP^/ VG (mi 
Using Eq. (|5T|) one finds 



m2) 



(60) 



h'^h, 



9ti^\p\^Lj^G^f3 



*i.^)+/.^ 



/3c 



(61) 



FIG. 1: The points are the numerical evaluation of fc{x, 2) + 
f^{x,z) (for three different values of 2 < 1) and the lines are 
plots of the corresponding fitted function given in Eq. (|63l) . 
Notice that the fitting function breaks down as we approach 
2 — >■ 1, which corresponds to 2lj — )■ /3c. 



Following the standard approach (see e.g., Ref. |l7J), one 
finds that the rate of energy loss from a system, in the 
far field limit, is given by 



dF 



20G 



|rPL,/i^^' . 



(62) 



This allows us to explicitly test the approximation 
against binary pulsar measurements, for which the en- 
ergy loss has been very well characterized |18| . 

For a quantitative fit to the data one would have to ex- 
tend this example to non-circular orbits and also account 
for tidal and other near field effects. However, such cal- 
culations rapidly become rather involved and even within 
this simple (and important) system one can derive some 
general consequences. 



The functions in Eq. (j58p are highly resonant at z = 
1, which corresponds to an orbital frequency uj — /3c/2, 
however they are readily calculated for both z < 1 and 
z > 1. In these regions the functions given in Eq. ([55)1 
can be evaluated numerically and fitted to an explicit 
functional forms. For w < /3c/ 2 this gives 



■ exp 



C 



f.imf- 



2lo" 



,/3|r|(l-^)"'V ■ cP, 



Ji /3|r| 



,(63) 



where C ~ 0.175 is approximately a constant except as 
2a; — )■ f3c. In Fig. [1] we illustrate some examples of this 
approximation and show that the approximation is good 
even for 2a;//3c sa 0.99. 
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FIG. 2: The points are the numerical evaluation of fc{x, z) + 
fs{x, z) (for three different values of z > 1) and the lines are 
plots of the corresponding fitted function given in Eq. (|64p . 
Notice that the fitting function remain a good approximation 
even as 2: — > 1, which corresponds to 2ui — >■ /3c. 



Whilst for a; > /3c/ 2 one finds 



/.|.H,|)]%[/.(^M,|^)- 



;2sinM/3|r|(/ 




{my 

where for the function / is approximately 



m 



(64) 




+2exp 



(65) 



Figure [D shows some examples of the quality of this fit, 
which, just as for Eq. (p5)) . is best away from a; — /3c/2, 
however remains very good even as one approaches this 
Hmit. It is important to note that Eq. (|64p does not 
approximate the General Relativity solution for /? — >■ cx), 
since, in this limit, all (finite) orbital frequencies satisfy 
Lo < /3c/2. 

Using the approximation given in Eq. (j63p . one can 
check that, for slow orbital frequencies, the expected 
result of General Relativity is indeed recovered in the 
/3 — >■ 00 limit. Specifically, we can expand Eq. (IS^ in 
the large distance (large |r|) and small orbital frequency 
(i.e., 2ui <^ c/3) limit, to find the first order corrections 
to the standard result of General Relativity, namely 



d£ 
'dt 



32G/i>- 



2„2, ,6 



5c5 

/3|r|(l-f^^ V Pc 



.(66) 



Thus, the f3 —> 00 (i.e., ap —> 0) limit reproduces the 
General Relativity result, as it should. Also note that 
any deviation from the standard result is suppressed by 
the distance to the source, at least for orbital frequencies 
small compared to /3c. Although the amplitude of the 
deviation from the standard result is small, there are two 
interesting features: firstly, the existence of a maximum 
frequency /3c and secondly, the fact that the deviation is 
oscillatory. 

The maximum frequency comes from the fact that, in 
addition to a natural speed c, this theory has a nat- 
ural length, given by /3^^ = -y/— aoG. This natural 
length scale comes from the first two moments of the 
test function used to define the spectral action, Eq. ^. 
Physically, one can think of this as the scale at which 
noncommutative effects become dominant. This is ex- 
tremely suggestive of an underlying maximum frequency, 
which would have rather significant consequences for the 
particle physics sector of the theory, in particular for 
renormalization. However, it must be remembered that 
here we are working with one simple system and even 
in this case Eqs. (|63p - (l64l) are numerical approximations. 
Nonetheless, the existence of a maximum frequency in 
this system allows NCG effects to significantly enhance 
the production of gravitational radiation. This is partic- 
ularly important given the suppression with l/|r|, that 
is present in Eq. (|66p. 

The presence of the Bessel function in Eq. ((66|) means 
that the amplitude of the deviation from the standard re- 
sult of General Relativity will oscillate both with chang- 
ing distances and changing frequencies. This allows for a 
myriad of possible observational signatures, such as dis- 
tinct beats of the observed energy loss of binary pulsars, 
correlated to their changing orbital frequency and the dis- 
tance to the binary. The fact that a similar phenomena 
occurs for the gravitational wave itself, Eq. ((5T|) , implies 
that there would also be a beat structure in direct de- 
tection observations. In the case of a binary pair the 
amplitude of the beat will be heavily suppressed com- 
pared to that of the carrier wave and thus it is likely to 
remain below observational sensitivity, except in systems 
with very large orbital frequency. 



V. CONCLUSIONS 

NonCommutative Geometry is a natural extension to 
our familiar notions of Riemannian geometry, that has 
the additional benefit of producing the action of all 
the Standard Model fields in addition to gravity terms, 
purely through geometrical considerations. Thus NCG 
treats both gravity and matter on an equal footing and 
provides us with concrete relationships between matter 
and gravitational couplings. The gravitational sector of 
(the asymptotic expansion of) this theory produces mod- 
ifications to General Relativity and in this paper we ex- 
plore the ramifications of these modifications on the for- 
mation and evolution of gravitational waves. 
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We have shown that the theory contains both massive 
and massless gravitons and that the requirement that the 
mass of these gravitons be positive fixes the sign of one 
of the couphngs in the theory (for a given choice of sign 
conventions). We also show that both these modes are 
sourced by the quadrupole moment of a system (just as 
in standard GR) and that the retarded Green's function 
is not restricted to the past light cone of the observer 
(unlike GR) , as one would expect for a system with mas- 
sive modes. We have explicitly calculated the energy loss 
for a circular binary system and compared the results to 
those of standard GR. We have demonstrated that the 
amplitude of these NCG modifications is suppressed by 
the distance between the observer and the source of the 
gravitational waves and hence will typically be small. 

Despite the extremely small amplitude of deviations 
from standard results, we have shown that NCG produces 
several distinctive features. Firstly, the amplitude of the 
energy lost by a binary pair can be higher or lower than 
the expected value, depending on the orbital period of 
the pair and the distance to the observer. This opens up 
the possibility that the observed energy loss from such a 
pair would be seen to oscillate as the binary moves with 
respect to the Earth. Whilst such effects are likely to 
be beyond current observational resolution, they allow 
for an unexpected beat phenomenon, which would be a 
concrete signature of NCG. 

In addition, we have shown that the amplitude of these 
effects is (approximately) proportional to (1 — 2uj/cfi)~'^ , 
where u is the orbital frequency of the binary. Thus, it 
would appear that the NCG corrections to the energy 
loss by the binary can become arbitrarily large as the 
frequency of the binary approach the critical frequency 
/3c. In such a regime, the weak field approach taken here 
would no longer be valid (and numerical approximations 
break down) , so one would not trust systems very close to 
this limit, however it is certainly true that astrophysical 
constraints on the parameters of the theory will be signifi- 
cantly improved for objects with a very rapidly changing 
quadrupole moment. A precise understanding of such 
systems is likely to require detailed knowledge of various 



astrophysical effects (radiation and particle production, 
tidal stripping etc.) as well as analytic solutions to the 
graviton field in the large frequency regime. 

Finally, the form of Eq. ([57]) suggests that similar be- 
haviour may be present in other systems, with periodic, 
or almost periodic, variations in the (mass) quadrupole 
moment. For laboratory systems, the gravitational radi- 
ation predicted by General Relativity is negligible, how- 
ever if the NCG enhancement were sufficiently large, this 
may no longer be true. Laboratory systems regularly 
have very large oscillation frequencies (e.g., lattice vibra- 
tions in solids can easily exceed lO^^Hz) which would ex- 
perience anomalous damping, if the system was produc- 
ing significant amounts of gravitational radiation. This 
opens up the (remote) possibility that the noncommuta- 
tive nature of space-time might be probed in the labora- 
tory. 

One can immediately use the results of this paper to 
examine circular binary systems, in order to constrain 
the value of /3 [10| . Similarly, one can include eccentric- 
ity which may result in more restrictive constraints on 
the theory. An alternative avenue would be to use the 
gravitational wave-forms given here to deduce the conse- 
quences for direct gravity wave searches (LIGO, VIRGO, 
LISA, etc). In particular, to extend these result to the 
large field regime and look for modifications to the chirp 
that develops at the end of in-spiral events. 
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FIG. 3: The retarded Green's function is defined by extend- 
ing the contour of the integral into the positive half of the 
imaginary plane around each of the four poles, to = ±]k] and 
a; = ±i/]k|2+/?2. 



Appendix 

In what follows, we shall explicitly detail the calcu- 
lations which were performed in order to determine the 
Green's functions, namely Eqs. (H71) . (|i5|) . 

Consider first the following fourth-order wave equation 

(D-l3^]DGnix,x') = inS^^^x - x') . (A-1) 

Lorentz symmetry of the background restricts Gr, to be 
solely a function of a; — a;'. Given the following inverse 
Fourier transforms: 



GTiix,x') = 



1 



'27r) 



d'*fcGR(fc)e''" 



ik-{x—x') 



and 



5{x 



{2nr 



j4^gjfe.(x-x')^ 



(A-2) 



(A-3) 



where k ■ z = 
must satisfy 

Gnik) = 



-OJZq 



k • z, the Fourier transform G^{k) 



in 



[{lo + ief - k2 - p^] [{oj + leY - k2] 



,(A-4) 



in order to solve Eq. (jA-ip . 



Upon performing the inverse Fourier transform to 
determine the coordinate expression for Gr(x,x'), the 
following pole prescription uniquely determines the re- 
tarded Green's function (the reader is referred to Fig. [3]). 
The contours are traversed above the poles in the com- 
plex u) plane: for i — i' > we close the contour in the 
lower half plane picking up the residue of the poles; for 
t — t' < Q we close the contour in the upper half plane, 
thus enclosing no poles. 



Note that GR,(fc) can be rewritten as 

1 



GR(fc) 



Att ■ 



1 



^ 13^ [(w -f iey - k2 - /32] 
1 

gi{k)-g2{k) , 



(A-5) 



which is in fact simply the Fourier transform of Eq. (^5 
We will analyze each term in Eq. (|A-5p separately. 

First we need to define the following quantities: 



T 
R 

k 



t~t' , 
r — r' , 

|k|2+/?2 



(A-6) 



In what follows, we will set the speed of light c — 1. The 
.gi(fc) term can be integrated as follows: 



gin{x -~ x') 



1 



(2^)4 

1 
4^ 



da;d^k.gi(fc) 
dwd^k 



„ — icjT^ik-R, 



[((w + ie) + \/k){{u} + ie) - Vk)] 



(A-7) 



where the -\-ie above is simply a mnemonic for the re- 
tarded Green's function pole prescription. The contour 
integral in the complex a; plane results in — 27r? y^ Res 
and we have 



gm{x - x') 



*e(T) 

47r2 



d^ke^''-^ 



^iT\fk _ -iT\fk 



Vk 



(A- 
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Upon performing the angular integral, one finds 

e(T)/3 



gm{x - x') = 



27r|R| 



[Z+-Z-] 



(A-9) 



where u € [uq, oo), uq = arcosh (T/t) and we have defined 



where 



Z± = 



fcdfc 



VF 



^/'g''/3(|R. 



\k±TVW+l) 



,-'i/3(|R|A:±TVP+T) 



(A-10) 



We will first focus on the solution which is interior to the 
light cone, i.e., T > |R|. In the first term in the inte- 
gral, Eq. (JA-9P above, we perform the following change 
of variables: 



\R\k + Ty/P + 1 = Tcoshw, 
Tfc+IRlVfc^ + l = rsinhu. 



dfc 



yPTT 



IR I T 
k = cosh V -\ sinh v , 

T T 

= dv , 



In the second term in the integral in Eq. (|A-9p above, 
we will perform a different change of variables given by 



|R|fc-TVA?2+T 

Tk- iriv'pTT 



—T cosh II , 
r sinh v , 

IRI 



T 



dfc 



VF 



■ cosh u H sinh ii . 

T T 

dv , 



where for these change of variables v € [— i;o,ck)). Note 
that, the variable v spans the interval [— tiQjO], while fc 
correspondingly spans the interval [0, |R|/t]. After these 
changes of variables we arrive at the following: 



gm{x - x') = 



e(r)/3 [ 

27r|R| \ 



e(T)/3 
27r|R| 

e(r)/3 

2t:t 



dv 



dv 



cosh V -\ smh v 

T T 

cosh u H smh w 

T r . 



2/3rcosht; i — z/^rcoshf 



2/3rcosht; _i_ — i/3rcoshij 



dw 



cosh w H smh v 

T T 



if^T cosYvv i —i j3r cos\vv 



dv 



i(3r coshv—v _|_ —i(3T coshv — v 



(A-11) 



r 



where the last equality follows by symmetry; it should 
be understood that the limits in the integrals above are 
such that the integrals are convergent. 

We note the following integral representations '22] of 
the Hankel functions of order a of the first and second 
kind, respectively: 



H^^\x) = 



.,— ia^ foo+ic 



m 
e 



dv I 



IX cosh v — av 



oo — ie 



1 —ix cosh v—av 



— oo+ze 



related to the Bessel functions of first and second kind 
via the relations: 



Hi^'>{x)^Mx) + iyc.{x), 



(A-12) 
(A-13) 



where Ja{x) and ya{x) are Bessel functions of the first 



and second kind respectively; i.e., the two linearly inde- 
pendent solutions to Bessel's equation: 



,(fF dF 
dx"^ dx 



{x^ ^ a^)F = . 



(A-14) 



In particular, we require the Hankel functions of order 
a — I, which are given as 



-1 />oo-t-2e 

h[^\x) = / dv 

TT J_oo-i£ 



^tx cosn v—v 



h['\x) 



oo—'ie 
oo — ie 



dv e 



-ix cosh v — v 



oO'\-ie 



Thus, one can express the integral for ^ir, Eq. (JA-Ill) . 
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in terms of these Hankel functions as 

e(T)/3 



giR{x-x) = 



27rr 



{7rHi'\M+^H['\(3T)} , 



= e(T)^^i^ for T > |R| . (A-15) 

T 

Note that, Ji(/3r) is the Bcsscl function of first kind of 
order 1 and we have used Eqs. (jA-12[) and (jA-13[) to arrive 
at the final expression, Eq. (jA-15|) . above. 

Looking now at the exterior of the hght cone, i.e., |R| > 
T, we make the following change in variables in the first 
term in the integral of Eq. (jA-9[) above: 



|R||k|+TV|kF+l ^ Csinhw, 
r|k| + |R|v/|k|2 + l = ecoshw, 

k 

d|k| 



^1 • V, ^ >, 
— smh V cosh v , 



, = dv . 

where as before v G [110,00), but now vq = arsinh (T/^) 
and we have defined ^ = a/|RP — T. 

In the second term in the integral of Eq. (|A-9p we per- 
form the following change of variables: 



T 

I 



^l^^ + l 

where v S [—vq, 00). 

Upon using these change of variables we arrive at 



|R||k|-rv/|kp + i 


= ^ sinh V , 


r|k| + |R|v/|k|2 + i 


= ^ cosh V , 


k 


|R.| . , 

— —— smh 


d|k| 


= dv , 



gmix - x') 



-'-i^h^ 



r2),(A-i6) 



which vanishes since we are explicitly considering the 
|R| > T region and thus the delta function vanishes. 

We still have yet to determine the singular part of the 
Green's function on the light cone |R| = T. To do this we 



where Gauss' theorem was used to arrive at the first term 
above and dS^ is the surface element of the boundary 
dV. Assuming Jy /i vanishes as the integration volume 
vanishes, we are left with 



(A-18) 
given by 



hm / (5i)^^dE^=47r. 

Now introduce the coordinates (w7,x, 0, 

t = w cos X , 

X = w sin X sin cos (j) , 

y — w sin x sin 9 sin </> , 

z — w sin X cos 9 , 

such that the line element ds^ = gaisdx^dx^ of the flat 
background takes the form 

ds = — cos 2xdw + 2w sin 2xdwdx 

+w'^ cos 2xdx^ + w^ sin^ xd^^ , (A-19) 



where dJl^ = d9^ + sin^ 



In these coordinates, the 



surface dV is given by constant w, and the Synge world 
function a is 



1 „ 

a = --W cos2x ; 



(A-20) 



notice that for time- like events — 2(t = t^ where r is as 
previously defined. 

The following quantities will be useful for what follows: 



'-9 



9 



XX 



3*2 

w sm X sm ' 
— cos2x , 
sin2x 

W 

cos2x 



where g is the determinant of the metric t/o,^ . The only 
nonzero component of dS^ is 



dS^ 



w^ sin xdxdi^ , 



(A-21) 



will repeat the formalism established in Ref. [2l[. Note where dfl ~ sin 



that, we very well could have determined the full Green's 
function giYi{x—x') and n ot .i ust the part on the light cone 
via the formalism in Ref. [2l|]. In fact, a full review of the 
formalism will serve as a useful check on the calculation 
of the smooth part of the Green's function determined 
above via Fourier transform. 

To begin, we first integrate the Green's function equa- 
tion, Eq. (|A-ip . over a space-time volume which contains 
the source event x' = 0, namely 



In these coordinates, the retarded Green's function is 
given by: 



giK = e(wcosx).9(cr) 



(A-22) 



where g{a) is an as yet undetermined, possibly distribu- 
tional, function. We will only need the following gradient 
of 51R (omitting the label r) 



av 



(51) 



ME 



^ /?' 



ffi = 47r , 



(A-17) 






(A-23) 
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A straight-forward calculation leads to 

(ffi)'"' = - S{w cos x) cos xgicr) 
+w<d{'w cos x)9'{'^) 7 



(A-24) 



where the prime on g{<7) denotes differentiation with re- 
spect to a. We then have 



dv 



(5i)^^dS, 



dV 



w sin xdxdfi { w9 (w cos x)ff' (c) 
(5(w;cosx)cosxff(CT)} 

4 /■' .,;„2,,J,,„/,„^ (^_25) 



47rw'' / sin x^Xg' (<^) , 

'0 



where the Heavyside function has restricted the limits 
of X integration such that cos x > and the delta term 
vanishes. 

Changing integration variable from x to cr in the inte- 
gral above we arrive at the following condition on g{cr): 



lim e 

£-i-0 



da 



H(^).V, 



1, 



(A-26) 



where 



« = 2«. 



-(7) 



'1 + f 



We now propose the following ansatz for g{(T): 

g{a) = V{a)Q{-a)+A5{a)+B5'{a) 

+C5"{a)+D5"'{a) + ... , (A-27) 

where V{a) is a smooth function of a and A,B,... 
are constants. Inserting the ansatz Eq. (|A-27[) into 
Eq. (|A-26P gives 



Imre^/ daS(^)y'(a)e(-a)-S(0)nO)-^; 



B, 



C_, 



(0) + ^S(0)-^S(0) + -SW(0) + 



r 



= 1 . (A-28) 



The first two terms on the left-hand side of Eq. (JA-28P 
above vanish, since V{<t) is assumed to be smooth. The 
limit exist as long as B = C = ... = and the condition 
Eq. (jA^ wiU be satisfied iff A == -1, since S(0) = 1. 
To fully determine the smooth part of g{u) one needs 
then only to solve the homogeneous equation (i.e., x ^ 



from which it is straight forward to verify that (the reader 
is referred to Ref. 21] for further details) 



V{a) = 



/3Ji(/3y^2^) 
\/^2^ 



(A-29) 



Returning to the original coordinates (t, r, 9, 0) we then 
have 



gm{x-x') = ^5{\t^)Q{T) 

r V2 

2\ 



Rp) e(r) 



+ 



2,5(r^)e(r) 



e (T - |R|) e(T) , (A-30) 



where we have used some standard properties of the Dirac 
delta function and the Heavyside step function. 

To determine g2Rix — x') we need to consider the in- 
tegral 



g2K{x - x') = — / dwd^k 



(A-31) 



„— iuT^Jk-R, 



[{uj + ie+\\i\){LO + ie~\k\)\ 



where the pole prescription is again that of the retarded 
Green's function. A straight-forward calculation leads to 



92k{x - x') == 



e(r) 

IRI 



{<5((T+|R|)-<5(T-|R|)}, 



-2e(r),5(7 



(A-32) 



Finally, the Green's function which satisfies the fourth- 
order wave equation Eq. (|42|) is given by 

1 



Gr^x-x') = -j^{gi^(x-x')- g2K{x-x')), 

-e(r-|R|)e(r). (A-aa) 



/32 
/3Ji(/3r), 



Pt 



Note that the Green's function above is subject to the 
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initial value condition: 

9Gr(x - x') 



dt i*=o ^ 



